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Abstract 

It is shown that at collider energies experimental distributions in the multi- 
plicity n of negatively charged particles in inelastic and non-single diffractive 
pp collisions are well parameterized by a sum of so-called Gupta-Sarma distri- 
butions having the Poisson distribution as a particular case. This extends the 
earlier description of the multiplicity distributions in hadron-hadron collisions 
at cm. energies below 65 GeV by the two parameter sum of Poissonians. Im- 
plications of the proposed parametrization for the LHC energy are discussed. 
PACS numbers: 13.85.Hd, 12.39.-x 
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I. INTRODUCTION 



Charged particle multiplicity distributions in inelastic hadron-hadron collisions at high 
energies are usually described by the negative binomial (NB), for an experimental review 
see Ref. [|l|. However the appearance of the shoulder structure observed for the first time 
by the UA5 Collaboration has led to the use of the weighted sum of two NB's with 
five free parameters M , where the first NB describes the contribution of soft events (events 
without minijets) and the second one describes the contribution of semihard events (events 
with minijets). The aim of this paper is to extend to the collider energies another phe- 
nomenological parametrization giving better agreement with lower energy data than NB 
both for p(p)p [|,f| and meson-proton collisions. In the Refs. @,[| the multiplicity 
distribution P n of negatively charged particles produced in inelastic p(p)p collisions at the 
center of mass energies \fs below 63 GeV have been fairly well described by a two param- 
eter sum of Poissonians. This approach is based on a simple minded quark-parton model 
in which quarks q interact pairwise independently of one another with the same conditional 
probability e and each qq interaction leads to the same multiplicity distribution in the final 
state. The probabilities for events with 0, 1, 2 or 3 qq interactions are equal respectively 
to (1 — e) 3 , 3e(l — e) 2 , 3e 2 (l — e) and e 3 and in terms of a probability generating function, 
p.g.f. (for a mathematical formalism see Refs. it leads to the relation 

G{z) = Y,PnZ n = H-e + e^z)) 3 , (1) 

where G(z) is the p.g.f. for the final distributon and f{z) is the p.g.f. for events with one 
parton-parton collision. The p.g.f. 's for events with two or three parton-parton collisions 
are simply convolutions tp 2 (z) and tp 3 {z). Good description of the experimental data has 
been obtained with ip(z) = exp (S(z — 1)), the p.g.f. for the Poisson distribution. This 
parametrization had the strong energy dependence of the parameter e |4j], more smooth 
energy dependence of the e was observed || when the Poissonian was replaced by the Poisson 
distribution truncated at zero multiplicity with the p.g.f. (p'(z) = (<p(z) — <^(0))/(l — </?(0)). 
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The need for the truncation at zero is explained for pp collisions by non-zero electric charge 
of the initial uu, ud and dd pairs leading to the reduced probability to have zero charged 
particles in the final state. This is not the case for pp collisions where uu and dd pairs (5 
out of 9 combinations) are neutral. 

The model [f§-|6|] is modified in the present paper in the two aspects: 
a) The Poisson distribution is replaced by the so-called Gupta-Sarma (GS) distribution 
g|n| with the p.d.f. 

= e -p f r.^rl ) • (2) 



.1 + r(l -z) 

The GS distribution having the Poissonian as a particular case at r = and known in 



mathematical statistics under the name of the Polya-Aeppli distribution |TT| has physical 
interpretation in the framework of different models |12|- |I^1 , see discussion in Refs. |15| , [16| . 
In the approach advocated by Biyajima et al. |j2j and Finkelstein [|13j the multiplicity 



distribution originates from the Poisson distribution of some clusters, each cluster obeys 
Furry- Yule (or truncated at zero Bose-Einstein) distribution, finally it leads to the p.g.f. (2). 
One can note that the same form (2) is valid when the p.d.f. for cluster decay distribution is 
a linear fraction (1 + A(1 — z))/(l+r(l — z)), usual for the theory of branching processes, 
Ref. [p/?]. In the Gupta-Sarma approach ||[10| the system after collision is viewed as one 
highly excited hadron emitting entity obeying the branching process with the probability per 
infinitesimal time At to produce k new particles proportional to X k At, where A is positive 
constant. The solution of the corresponding evolution equation for initial condition with 
zero particles leads to the p.g.f. of the form (2). In the Chau- Huang approach |TJ| the GS 



distribution is obtained from the statistical Ising model. 

b) As suggested in the Ref. H events with zero parton-parton collisions can represent 
the diffractive-like processes, with the fraction of diffractive-like events given by the (1 — e) 3 . 
It has been established that the multiplicity distribution for diffractive system with effective 
mass M looks like the multiplicity distribution in pp collisions at the cm. energy ^Js = 



M, see Ref. [18| for a review on diffraction. In this paper we approximate the diffractive 



contribution by the form (1) with <p{z) equal to the p.g.f. for the Poissonian under the crude 
assumption that the integrated over M distribution is similar to the distribution at some 
effective mass M and the final p.g.f. is given by 

G(z) = (1 - e) 3 ip d + 3e(l - ef g(z) + 3e 2 (l - e) g 2 (z) + e 3 g 3 (z) (3) 

with 

<p d (z) = {l-e + e exp( S d (z - l))) 3 . (4) 

More careful description of the diffract ive- like events with integration over M is given in 
Refs. 



In the section II the main characteristics of the Gupta-Sarma distribution are summa- 
rized. In the section III we present results of fits to the available pp data and to the pp 
data obtained at the SppS collider. In the section IV the discussion and the conclusions are 
given. 



II. CHARACTERISTICS OF THE GUPTA-SARMA DISTRIBUTION. 

The mean multiplicity < n > and its dispersion D = (< n 2 > — < n > 2 ) x / 2 are 
easily obtained from the p.d.g. (2) using formulae < n >= dg/dz \ z= \ and D 2 =< n > 
+d 2 In g(z)/dz 2 | 2=1 

<n>=S, D 2 =< n > (1 + 2r). (5) 

To obtain expressions for probabilites g n one can use the method proposed by Finkelstein 
||13|| . The p.g.f. is expressed as a sum of powers (z/(l + r (1 — z))) k , where the denominator 
represents the well known NB. Then, the contribution to g n with n^O from the fc-th term 
is equal to the NB probability to have n — k particles and finally it gives 

with 
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and 

g = g(0) = exp (-^-) . (8) 



One can note |T2|JT5| that the GS distribution is a particular case of the partially coherent 



laser distribution PCLD (see review of the PCLD in 0]). Indeed the p.g.f. for the PCLD is 
the product of the p.d.f. (2) and (1 + r (1 — z))~ k , i.e. the convolution of the GS and NB 
distributions. It gives the expressions for g n in terms of the Laguerre polynomials [^1[ 



9n = (iTtY exp (i^) l "(t(iT7)) • (9) 



Using iteration relations for the Laguerre polinomials [^TJ one can obtain next iteration 
relations for g n 

(n + 1) g n+1 = (a + 2n b) g n - (n - 1) b 2 g n - X (10) 

at n > 1 and 

9i = ag - (11) 
These iteration relations can be useful for calculations at large n values. 

III. RESULTS OF FITS 

Both for pp and pp data we calculate the number of negatively charged particles as 
n = (n c h — 2)/2, i.e. we count the number of produced pairs of charged particles. As 
mentioned above, we truncate the distribution for parton-parton collision in the fits to the 
pp data and do not truncate it for pp data at the collider energies. 

The pp data used [22-39] are the same as in Refs. ]4]|| with additional measurement 
from Ref. [3~I|. In the Table 1 the results of fits to the distribution with the p.g.f. (3) are 



given for the case r = 0. The agreement with experimental data is good, this is expected 



since even the two-parametr parametrization with the diffractive contribution concentrated 
at zero multiplicity was successful [|]||. One can note also that the mean multiplicity for 
diffractive contribution, proportional to Sd increases slowly with energy. 

Non-single diffractive (NSD) multiplicity distributions, measured by the UA5 Collabora- 
tion [Q,fH]] have been parameterized by the distribution (3) without diffractive component, 
i.e. by the distribution with three other "parton-parton collision" components normalized 
by the factor (f — (f — e) 3 ) -1 . Inelastic multiplicity distribution at ^/s = 546 GeV f4(J has 
been parameterized by the full distribution (3). The results of the fits are given in the table 2 
both for e = 0.456 and for free e. The fixed value of the e was chosen on the assumption that 
the fraction of diffractive like events is equal to 16%, the fraction of the single diffractive 



events measured by the CDF Collaboration at 546 GeV |4I[] . The fraction measured by the 
UA5 Collaboration is equal to If % f40[l , corresponding conditional probability e is equal to 
0.52. 

The results of the fits with free e are illustrated in Figs, f-3 respectively for inelastic 
and NSD data at 546 GeV |0J and for NSD data at 900 GeV §. The quality of the fits 



is quite qood, the fluctuations in the parameter e are explained by the our crude treatment 
of the diffractive component and possible bias in the experimental data. We have ignored 
also the nonnegligible contribution of the double diffraction processes in the NSD data. The 
influence of the high multiplicity tail on the fit parameters has been observed also, the fit in 
the region n c h < 80 of the inelastic data at 546 GeV gives more reasonable value of Sd near 
3, significantly smaller than the values Sd in the Table 2. 

IV. DISCUSSION AND CONCLUSIONS 

The possibility that multiplicity distributions at high energies can split into several struc- 
tures has been predicted more the 25 years ago E3] . For example Nielsen and Olesen made 



the statement jKJ: "if we go to high enough energy one should see a separation of the 
multiplicity spectrum in a series of equidistant peaks at n ~ n\, 2ni, 3ni, Kaidalov 
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and Ter-Martirosyan in the framework of the quark-gluon string model have predicted three 
peaks in the multiplicity distribution at y/s = 100 TeV f44j . These Regge-type models 
P^ , |44|| in principle predict more than three peaks (structures) in contrast to our approach 
with maximum three nondiffr active structures. The prediction for NSD multiplicity distri- 
bution at LHC energy y/s = 14 TeV, calculated with parameters e = 0.456, r = 0.8 and 
5* fixed by the expected mean multiplicity < n c h >= 67.2 is given in the Fig. 4. It seems 
to be intermediate between the drastic predictions of the Regge-type models and more flat 
predictions based on the parametrization by the weighted sum of the two NB's [[45| , f46|l . 

In conclusion, the multiplicity distributions at collider energies have been fairly well pa- 
rameterized by the sum of Poisson-like GS distributions, with one GS distribution describing 
the multiplicity distribution for one "parton-parton" collision. An attempt has been made 
to connect the fractions of events with 1, 2 or 3 "parton-parton" collisions with the fraction 
of diffractive-like events in a framework of the simple minded parton model. 
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TABLES 

TABLE I. Results of the fits to the negatively charged particle multiplicity distributions in 
inelastic pp collisions [22-39]. 



Ref. 


v 7 * (GeV) 


e 


S 


s d 


X 2 /NDF 


@ 


6.84 


0.336±0.002 


0.203±0.006 


0.032±0.013 


21.1/5 


@ 


7.87 


0.381±0.004 


0.325±0.011 





3.9/5 


@ 


9.78 


0.441±0.015 


0.521±0.040 





7.5/6 


m 


10.69 


0.470±0.014 


0.500±0.044 





2.8/5 


m 


11.46 


0.470 ±0.008 


0.671±0.019 


0.013±0.040 


18.6/7 


m 


13.76 


0.472±0.009 


0.940±0.028 





19.1/7 


m 


13.90 


0.483±0.010 


0.855±0.032 


0.003±0.041 


7.3/7 


m 


16.66 


0.491±0.006 


1.170±0.021 


0.012±0.021 


18.3/10 


u 


18.17 


0.523±0.031 


1.224±0.094 





2.4/7 


a 


19.42 


0.551±0.016 


1.197±0.054 


0.074±0.074 


9.0/8 


n 


19.66 


0.538±0.011 


1.243±0.034 


0.070 ±0.071 


7.3/10 


n 


21.7 


0.512±0.011 


1.460±0.040 


0.135±0.063 


14.3/11 


a 


23.76 


0.563±0.012 


1.486±0.039 


0.112±0.066 


10.4/11 




23.88 


0.561±0.019 


1.601±0.074 


0.246±0.148 


11.2/10 




26.0 


0.577±0.010 


1.638±0.034 


0.102±0.054 


8.4/10 




27.6 


0.542±0.013 


1.720±0.059 


0.096±0.057 


17.8/13 




27.6 


0.555 ±0.016 


I p* r* f i r\ r\'~7'~7 

1.565±0.077 


0.027±0.068 


4.3/9 




30.4 


0.528±0.021 


2.000±0.077 


0.300±0.157 


3.1/14 


@ 


38.8 


0.576±0.008 


2.059±0.031 


0.273±0.085 


8.3/13 




44.5 


0.538±0.021 


2.441±0.086 


0.438±0.167 


4.9/16 




52.6 


0.549±0.016 


2.647±0.068 


0.381±0.142 


12.2/18 




62.2 


0.552±0.017 


2.887±0.070 


0.426±0.168 


16.7/17 
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TABLE II. Results of the fits to the negatively charged particle multiplicity distributions in pp 
interactions [2,40]. 



events 


(GeV) 


e 


S 


s d 


r 


X 2 /NDF 


NSD 


200 


0.456 


5.965±0.076 




0.170±0.056 


19.7/29 






0.264±0.038 


7.441±0.304 




0.297±0.084 


9.2/28 


NSD 


546 


0.456 


8.429±0.053 




0.468±0.030 


61.3/45 






0.352±0.018 


9.453±0.187 




0.558±0.038 


32.6/44 


inel. 


546 


0.456 


7.410±0.091 


11.071±0.458 


0.832±0.075 


39.4/44 






0.536±0.027 


6.743±0.067 


11.811±0.326 


0.794±0.067 


29.9/43 


NSD 


900 


0.456 


10.400±0.100 




0.703±0.065 


77.0/52 






0.304±0.029 


12.213±0.394 




0.823±0.095 


20.9/51 
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FIG. 1. Multiplicity distribution for inelastic antiproton-proton collisions at the cm. energy 
546 GeV [40] (squares) compared with results of the fit (full dots). 
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FIG. 2. Multiplicity distribution for NSD collisions at the cm. energy 546 GeV [40] (squares) 
compared with results of the fit (full dots). 
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FIG. 3. Charge particle multiplicity distribution for NSD collisions at the cm. energy 
900 GeV [2] (squares) compared with results of the fit (full dots). 
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FIG. 4. Prediction of the NSD charged particle multiplicity distribution at the LHC energ 
14 TeV. 
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